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Density  functional  theory  has  been  an  important  method 
for  studying  the  ground  states  for  systems  involving 
interacting  electrons  such  as  atoms,  molecules,  and 
solids.  This  paper  helps  to  extend  this  theory  to  excited 
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expressions  for  these  quantities,  valid  for  systems  of  slowly  varying 
density.  Even  for  such  systems,  however,  there  are  essential  non-local 
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effects.  Nevertheless  both  E  and  v  can  be  calculated  In  terms  of 
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and  allows  calculation  of  excited  state  energies  and  densities. 
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ABSTRACT 

The  starting  point  of  this  paper  is  a  recent  extension  by  A.K.  Theophi- 
lou  of  the  Hohenberg-Kohn-Sham  (HKS)  density  functional  theory  to  ensembles 
of  systems  consisting  of  the  M  lowest  eigenstates,  equally  weighted.  As  in 
the  HKS  theory  the  key  quantities  are  the  exchange  correlation  energy, 
E^^Cn(rO]  and  potential  (r;Cn(r')]).  The  present  paper  provides  expres¬ 
sions  for  these  quantities,  valid  for  systems  of  slowly  varying  density. 

Even  for  such  systems,  however,  there  are  essential  non-local  effects, 
nevertheless  both  eJ^  and  v  ^  can  be  calculated  in  terms  of  quantities 
characteristic  of  appropriate  uniform  thermal  ensembles.  This  theory  is 
the  analog  of  the  ground  state  local  density  approximation  and  allows  cal- 
culation  of  excited  state  energies  and  densities.  ^ccesion  For _ Y 
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The  ground  state  density  functional  theory  of  Hohenberg,  Kchn 

and  Sham  (HKS)  '  '  was  formally  extended  by  Theophilou"'  to  the 

mean  energy,  E^,  and  mean  density,  n{r) ,  of  an  ensemble  consisting 

of  the  lowest  M  states  equally  weighted  by  the  factor  1/M.  (We 

shall  call  such  an  ensemble  an  equi-ensemble . )  This  theory  allows, 

in  principle,  the  calculation  of  individual  excited  state  energies 

and  densities,  E  and  n  (r) .  In  analogy  with  the  HKS  theory,  the 

m  m 

essential  required  quantities  are  the  exchange -correlation  energy 

and  potential,  E^  [n(r’)  ]  and  {r,[n(r')  ])  ,  both  functionals  of 

xo  «cc 

the  ensemble  density  n(r'). 

The  practical  usefulness  of  the  HKS  ground  state  theory  has 
been  largely  due  to  the  simplicity  and  surprising  accuracy  of  the 
so-called  local  density  approximation  (LDA) , 


E^^[n(r')  ]  =  i'e^^(n(r'))dr' 


clnv®xc'^*  n*n;r) 


)  LDA 


per  unit  volume 

where  is  the  exchange-correlation  energy^of  a  uniform 

electron  gas  of  density  n  in  its  ground  state.  The  present 
paper  reports  briefly  a  generalization  of  the  LDA  to  the  equi- 
ensemble.  A  complete  manuscript  has  been  submitted  to  Physical 


Review  A. 

We  remark  first  of  all  that  expressions  for  E^  and  v^ 

^  xc  xc 

analogous  to  (1)  and  (2)  do  not  exist.  The  reason  is  that, 

for  oi"en  >!,  the  contribution  of  a  volume  element  dr  to  E^ 

xc 
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aip^Ji-TOoogly  not  only  ott  th.  locil  on  the 

density,  n(r'),  everywhere  else.  Nevertheless,  a  quasi-LDA 

possible  for  the  equi-ensemble .  This  approximation  is  ohe  more 

accurate  the  larger  the  number  of  particles,  N,  ana  the  smoother 

the  density  n(r) .  The  derivation  proceeds  via  the  thermodynamic 

equivalence  of  such  a  large  and  smooth  equi-ensemble  with  a 

canonical  ensemble  of  appropriate  temperature  r.  Temperature 

ensembles  were  first  discussed  by  Mermin,"^^  and  local  density 

approximations  for  E®  and  v®  ,  quite  analogous  to  (1)  ana  (2) 

2) 

^  exist.  However  the  temperature  0,  of  the  canonical  ensemble 
equivalent  to  the  equi-ensemble  of  M  states  and  average  density 
n(r'),  depends  both  on  M  and  on  n(r')  for  all  r' .  Furthermore  the 
temperature,  6  ,  of  the  non-interacting  (KS)  canonical  ensemble 

which  is  thermodynamically  equivalent  to  the  non-interacting  (KS) 
ensemble  of  M  states  is  different  from  «. 

Our  final  results  are  the  following? 

Consider  a  system  of  N  electrons  in  a  given  external  potentia 
v(r).  The  objective  is  to  calculate  the  average  density,  n(r) ,  an 

w  C  \ 

average  energy,  E  ,  of  the  lowest  M  eigenstates  ' . 

1.  One  requires  the  following  thermodynamic  functions, for 

homogeneous  interacting  and  non-interacting  electron  gases, 
of  the  density  n  and  temperature  6'.  (The  subscript  s  denotes 
non-interacting  and  the  subscript  1  differentiation  with 
respect  to  temperature.) 

p  '  p  ' 

a)  The  entropies  per  unit  volume  a"  (n)  and  (n) ,  and 

0  ’  ^  * 

their  temperature  derivatives,  (n)  and  j_  , (n) . 

b)  The  exchange-correlation  plus  )cinetic  energy  per  unit 

pi  T  I 

volume  e  (n) ,  and  its  temperature  derivative  e^  (n) . 
c,  The  kinetic  energy  per  unit  volume  of  a  non-interacting 

,1  •  i  « 

svstem,  t„  (n) ,  and  its  temperature  derivative  t'  ,fn). 
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Begin  with  an  initial  approximation  to  n(r) .  Determine  the 

corresponding  interacting  temperature  9  and  non-interacting 

(KS)  temperature  9  by  solving  respectively  the  implicit 

s 

equations 


k  log  M  «  jTa®  (n(r)  )  dr;  k  log  M  •  (n  (D  )  dr. 
3 .  Construct  the  effective  one  particle  potential 


^eff^^’  *  j  fr^T  +  v”^(r?[n(r')]). 


where  v”  is  given  by  the  foLlowina  expression: 
xc 
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-  t  ®(n)j  /  X  +  4^/  \  I  e®  (n(.r')  )dr' 

s  ^  /n“n(r)  3n(r)  J  1 


30  r  9 

-  s  t  S  (n(r')  )dr’ 
3n(r)  /  s,l 


where 

30 _ 

3n(r) 


UA  /  :  a®(n(r’))dr' 

I  3n  ;n=*n(r)/  1 
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3n  (rl 


3a  ®(n)\ 

S  ) 


I -  ’  ,  (n(r')  )dr' 

3n  n=n (r)/  /  s,  1 
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4.  Solve  the  KS  single-particle  equations 


5. 


7. 


(-  +  Vg5f(r)  -  e^.)  (D'J'Cr)  •  0  .  (10) 

Construct  the  M  lowest  non-interacting  N-particle  wave- functions 
and  calculate  their  average  density  n'(r)^\ 

s  e«n 

If  n’li")  2  n(r),  then  the  original  n(r)  was  self-consistent.  If  not 
repeat  steps  2-5,  starting  with  a  different  initial  density  until 
sel f-consistency  is  achieved. 


M  -  - 

Now  determine  the  average  energy,.  £  ,  of  the  equi-enser.oj.e 

as  follows.  Let  ^(m=l,...M)  be  the  energies  of  the  M  . 

s ,  m 

lowest  KS  states.  Then 


M 

S*  =  Av 


ar 


r.  a 

( le" (n{r) )  - 


n  ir ; 


■ar 


(9) 


By  successive  calculations  for  increasing  M,  starting  with 

M  =  1,  the  excited  state  energies  E  and  densities  n^(r)  (averaged 

in  m 

over  multiplets)  can  be  obtained.  The  quantitative  accuracy  of 
this  quasi-LDA  for  systems  in  which  N  is  not  very  large  and/or 
r.  r'  is  not  slowly  varying  remains  to  be  tested. 
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